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Abstract We present the analysis of Spitzer/IRAC and NIR imaging observation of 
AFGL 5157, an active star forming region. In the IRAC images, this region shows strong 
PAH emissions in channel 4 and H 2 emissions in channel 2. Many of the H 2 features 
are aligned to form jet-like structures. Three bipolar jets in the NH3 core region and a 
couple of jets northwest of the core have been identified. We identify the possible driv- 
ing agents of the bipolar jets and show them to be very young. An embedded cluster has 
been detected in the NH 3 core; many members in the cluster show their SEDs increasing 
from JHK bands toward longer wavelength, indicative of their early evolutionary stages. 
Millimeter and sub-millimeter continuum emissions are found to coincide spatially with 
these presumable Class 0/1 sources, in the NH3 core and the NW subregion. The exis- 
tence of H 2 bipolar jets and very young stellar objects suggests that star formation is still 
going on at present epoch in these subregions. Combining the information from previous 
studies, we propose a sequential star formation scenario in the whole AFGL 5157 region. 

Key words: ISM:individual (AFGL 5157) — ISM: jets and outflows — ISM: lines and 
bands — stars: formation 



1 INTRODUCTION 



AFGL 5157 is an active star forming region that has attracted lots of attentions recently. Near infrared 
(NIR) observations show that a cluster of tens of young stars (maybe pre-main sequence stars) with a 
total mass of ~ 100 M are just emerging out of their natal cloud ( Chen et al.|1999 Kumar et al. 2006 



hereafter NIR cluster). An IRAS point source (05345+3157) is associated with the brightest NIR source 
IRS 1 in the NIR cluster. The far infrared luminosity of the IRAS source is estimated ~ 5.5 x 10 3 L Q at 
a distance 1.8 kpc ( |Snell et al.|1988|l, and this value could be scaled up if the recent kinematical distance 
2.1 kpc is adopted {Molinari et al.|2008| ). 

More interesting is the region ~ 1.5' to the northeast of the NIR cluster. NH3 emission lines were 



first reported by |Torrelles et al. ( 1992) and we refer this regio n as NH3 core her eafter. A number o f 
molecules tracing high density gas, such as HCN (pirogov 1999 1, N 2 H + and N 2 D + ( Fontani et al. 2008 ), 



have been detected recently. Millimeter continuum peaks have been detected in the core region extending 



to the south (|Klein et al.|2005|l. A molecular outflow was reported in the east- west direction (Snell et al. 



1 1988[ , Zhang et al. 2005), which has been resolved by the Submillimeter Array (SMA) observations into 
several outflows with a complex morphology (Fonta nTet al.|2009] l. A number of shocked H 2 emission 
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knots have been detected ( |Chen et al.|2003| l. Some of them are associated with the HH objects ( jTorrelles] 
|et al.| 1992] > and the outflows. These shocked H 2 emission knots, although are widespread and oriented 
towards different directions, show strong sign that they are associated with the NH3 core, suggesting 
the presence of a protostellar cluster with on going star forming inside. However, few NIR sources 
have been found in the region, suggesting that this proto-cluster is still deeply embedded in its nursery. 
We therefore carried out an observation using the Spitzer Space Telescope (SST) to investigate the 
embedded population and their properties. 



The observation reveals a large amount of deeply embedded sources that are not detected by Chen 
|et al.| ( |1999| l in the NH3 core, and a number of H2 emission knots being associated with the deeply 
embedded cluster. In this paper we present the highlights of the study. 

2 OBSERVATION AND DATA REDUCTION 

The observation was carried out on 2009 April 22, with IRAC mounted on the SST in imaging mode. 
Data in four channels (with central wavelengths 3.6, 4.5, 5.8 and 8.0 /im, hereafter chl, ch2, ch3 and 
ch4, respectively) were obtained. For each channel, a total of 288 frames were taken, each having an 
effective exposure time of 10.4 sec. After rejecting unusable frames (~ 30 frames in each channel), 
corrected basic calibration data sets were used to make the mosaic images using the software Mopex. 
Standard parameter inputs were used except for the fiducial image pixel scale, which was set to 0.72" per 
pixel, in order to obtain higher photometric accuracy. The final mosaic images have ~ 40 min integration 
time resulting in a 5er detection limit of ~ 0.1 mJy in all channels. 

The aperture photometry was carried out using DAOFIND and PHOT tasks in the IRAF package. 
Since the core region is rather densely populated with point sources, we adopted 4 pixels (2.88") as 
the aperture radius to do the aperture photometry. The fluxes were then corrected to 10 pixels (7.2") in 
chl and ch2, and to 15 pixels (10.8") in ch3 and ch4 because of larger expected point-spread-function 
in these two channels. The correction factors were obtained by comparing the fluxes of stars in the 
reference fields, which were observed simultaneously with the targets in each channel, between the two 
apertures (i.e., 2.88" and 7.2" in chl and ch2, and 2.88" and 10.8" in ch3 and ch4). The final factors 
are 1.21, 1.21, 1.38 and 1.42 for the four channels, respectively, with uncertainties better than 10%. The 
corrected photometric data were then combined with the 2MASS point source catalogue to obtain the 
SEDs of young stellar objects (YSOs). 

3 RESULTS AND DISCUSSIONS 

Fig.[T]presents pseudo-color image composed of chl, ch3, ch4 (upper panel) and chl, ch2, ch3 (lower 
panel) frames, with the blue color corresponding to the shortest wavelength and the red to the longest. 
In the upper panel, strong red extended emission is seen in the whole region, roughly extending in the 
southwest-northeast, indicative of the existence of PAH molecules. A shell-like structure surrounding 
the NIR cluster ( |Chen et al. 1999) is clearly seen. Note that the feature extending in the northern- 



southern direction might be artifacts due to the strong emission associated with the shell structure. We 
indicate them with arrows in the upper panel. To the northeast of the shell structure, in the NH3 core 
region, diffuse emission features are highly structured, similar to the M17 region (|Povich et aL| 2007 ) 
and other star forming regions (e.g. AFGL 437 Ku mar Dewangan & Anandarao||2010[ l. However, by 
careful inspection we have not found a sharp edge of the PAH emission as in the case of Ml 7 ( |Povich| 
et al. 2007). This suggests that the AFGL 5157 region lacks of very high-mass stars (e.g., O stars), thus 
lacks EUV photons to destruct the PAH molecules. 

In the lower panel, extended emissions are also obvious, similar to the upper panel. In the NH3 core 



region (Torrelles et al. 1992 ), a number of red point sources are detected. Many of them are not detected 
in the Ks-band and shorter wavelengths by the 2MASS and |Chen et al.| (|T999), suggesting they are still 
deeply embedded. The most prominent features are the compact and extended emissions in ch2, which 
appear green in the figure, generally close to the NH3 core region. These objects, commonly referred to 
as extended or compact green objects, are likely shocked H2 emission (Cyganowsk i et al.|200 8 ) arising 
from the star forming activities. 
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Fig. 1 Pseudo-color image composed from chl, ch3, ch4 (upper panel) and chl, ch2, ch3 
(lower panel) images, with the red, green and blue colors corresponding the decreasing wave- 
lengths. The R.A. and Dec. are given in J2000.0. 



3.1 The H 2 jets and knots 

To hightlight the H2 emission features, we present in Fig. [2^ the ch2/chl image, which is shown in a 
smaller field of view than that of Fig.[T] Although ch2/chl image can protrude the H2 emission features, 
some point sources that show stronger continuum emissions in ch2 than in chl cannot be removed 
completely. These point-like sources are circled in the figure. For comparison, we overlay line wing 
emissions of CO J=2-l observed with the SMA. The red and blue contours are integrated in tensities 
over velocity intervals (-9.2, 6.4) km s _1 and (-46.4, -26.0) km s _1 , similar to that presented by |Fontani 
[eTaLl ( |2009l l. 

First look at Fig. 2a reveals a large number of H2 emission knots in the NH3 core region and in the 
region ~ 2' west (hereafter AFGL 5157 NW). In the NH3 core region, the H2 emission features are very 
complex. The fact that tens of embedded YSOs exist in the region further complicate the picture. To 
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Fig. 2 a. ch2/chl image in a small field than that of Fig. [T] Some point sources having stronger 
emission in ch2 than in chl are circled for readers' reference. The contours are integrated 
intensity of the redshifted (red, -9.2 - 6.4 km s _1 ) and blueshifte d (blue, -46.4 - -26 .0 km 



s 1 ) components of CO J=2-l line emission observed with SMA ( jFontani et al.|2009 i. The 



contour starts from 4er (rms) with a step of 4<r (red: ~ 7.2 Jy beam km s ; blue: ~ 6.0 
Jy beam -1 km s _1 ). The dashed square represents the field of (b). b. Color image composed 
from ch2 (red), chl(green) and H2 narrow band (blue) images. The H2 narrow band (2. 12/im) 
image is adopted from Chen et al. ( 2003) >. The red contours are the 3 mm continuum observed 
by|Fontani et al. (2009 1. The green arrows and lines denote the rough direction of the bipolar 



jets in the area. The red pluses indicate the position of YSOs mentioned in the text. 



identify these H2 knots more concretely, we present in Fig. [2}} a color image composed from ch2 (red), 
chl(green) and H2 (2.12 fj,m) narrow band (blue) images. The image of H2 narrow band is adopted 
from |Chen et al. (2003 1. In this figure, the H 2 emission features are clearly seen in purple. By careful 
inspections, we could identify three bipolar structures and also a jet-like feature. Presumably they are 
jets arising from protostellar objects in the NH3 core. 

In the most prominent jet labelled as Jl, two bow shocks are seen in the two terminations, with their 
openings opposite to each other. In the midway of Jl, a point source is detected in all four channels 
as well as in the JHKs bands by 2MASS, being located at Yl: 05' l 37 m 52.02 s +32°00'04.1" (J2000). 
The red- and blue-shifted CO components (Rl, Bl, Fig. [2^) are surely associated with Jl. Rl has two 
peaks: one coincides the northeast H2 bow shock (Rla), the other is near Yl (Rib). Its blue-shifted 
counterpart Bl splits into several parts, denoted as Bla and Bib. Inspection of the CO channel maps 
indicates that Rib, Bla and Bib are high velocity components (~ ±25 km s _1 with respect to the 
system velocity ~ —19 km s _1 ), while Rla has a lower velocity (~ +15 km s -1 ). Remarkably, Bla, 
which is also split into two parts, is located at the wake of the western H2 bow shock of Jl, suggesting 
that the CO outflow is entrained by Jl ( [Canto & Raga| 1991 [[Riga & Cabritp993) . WISE (Wide-field 
Infrared Survey Explorer) archive image in the field shows a very strong mid-infrared point source at 
the Yl position, suggesting that this object is very young. Using the photometric data fro m the J-band 
to ch4 and WISE point source photometr)^] we fit the SED to the YSO models given by Robitaille et 



[aT]([2b07 1^] The best-ten fits give the mass of 2-5 M© and the age of 1.5-3.5 x 10 5 yrs. Fig! [3] shows the 
result of the best fit. We notice the strong excess emissions toward wavelengths longer than 10 /im. This 
is likely caused by contamination of other sources or interstellar dust emission since the FWHM of Yl 
in the fourth band (22 /im) of WISE is as large as 20". We also estimate the dynamical age by assuming 
the jet velocity of 25 km s _1 (the velocity of Bla) and th e inclination angle o f ~ 60°, resulting in an 
age of 2.2xl0 4 yrs at the assumed distance of 2.1 kpc ( Molinari et al.||2008 i. We note however, the 
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inclination angle is rather arbitrary, we cannot estimate the error arising from the assumption. The only 
meaning of the estimation of the dynamical age is that the driving source is quite young, consistent to 
the result from SED fitting. 

The second bipolar jet (J2) is inferred from the bow-shock at the northwest terminal and some faint 
H 2 emission knots in the pathway. At the southeast terminal, the morphology is irregular, probably 
due to the contamination of point sources, but we still notice a slight elongation in the jet direction. 
In the middle of J2, there are two point sources, located at Y2: 05 ,l 37 m 49.83 s +31°59'48.2" and Y3: 
05 h 37 m 50.02 s +31°59'52.4", respectively. The positional alignment with respect to the jet cannot sug- 
gest which one is the driving source. We then make the SED fittings to these two objects. The best fits 
suggest that both Y2 and Y3 are massive YSOs of ~ 8M , while Y2 is much younger (~ 1.6x 10 4 yrs) 
than Y3 (~ 3.8x 10 6 yrs). Fig.[3](right panel) shows the best fit for Y2. Since the point source detected 
by WISE is positionally associated with both Y2 and Y3, we assign the fluxes in four bands to Y2, while 
keeping in mind that they could be the joint contribution of both sources. This result would suggest Y2 
is more likely the driving source of J2. We have not found any high velocity CO component associated 
with J2, probably because it is near the edge of the SMA field of view. 

The third bipolar jet (J3) has a bow shock structure in the southeast and some faint tails in the 
northwest (Che n et al.|2003) >. The red- and blue-shifted CO emissions (R2, B2) are also detected along 
it (Fig. [2^). B2 is spatially coincident with the northwest H2 knot and R2 is a little closer in projection 
than the southeastern bow shock. The CO channel maps suggest R2 and B2 have veloc i ties ~ ± 15 
km s _1 . Between the two ends, no point source is detected in the infrared. |Fontani et al. (2008 2009 1 
observed this region with the SMA and PdBI. In their 1 mm and 3 mm continuum images, there is a faint 
peak northeast of CI, which is also associated with southern part of the northern N2H + condensation. 
As can be seen in Fig. |2j), the faint peak is in the middle of J3. Probably this mm source is the driving 
source of J3. Assuming a shock velocity of 15 km s _1 (the velocity of R2 and B2) and an inclination 
angle of ~ 60°, the dynamical age of this jet is estimated to be ^1.8x 10 4 yrs. 

Apart from the three bipolar jets, another jet-like knot (J4) is located between J2 and J3, elongated 
in the direction similar to that of J2. Faint H2 emission features have been detected in ch2 as well as 
in the NIR ( |Chen et aL"|20 03 ). However, we have not found any point source along the path. It is not 
possible to conclude where this jet is originated at present. A possible interpretation is that it is a parsec- 
scale jet from the NW region. In addition, there are a number of faint H2 emission features. They are 
distributed so disorderedly that we cannot tell from where they originate presently. Notably, there is 
likely another bipolar CO outflow with velocity ~ ±15 km s _1 (denoted as R3, B3) in the region, while 
R3 is contaminated by Rl and B3 is positionally connected to B2. 

In the NW region, two jet-like structures are found. One is stronger, roughly in the east-west direc- 
tion. Some knots have also been detected by |Chen et al.| ( 2003] l. The other fainter jet feature, oriented 
in the southeast-northwest, is not detected by jChen et al. (2003). These two jets have a number of com- 
ponents that are too faint to determine their exact direction, so we only place two dotted lines along 
these components for illustration. Nevertheless, they are likely originated from the region where they 
intersect, and where a number of point sources with red SEDs are detected. These red point sources are 
reasonable candidates of the driving sources of the two jets. High-resolution observations of molecular 
outflows to this region may help to understand the star forming activity there. 

In conclusion, the existence of jets/outflows in these subregions signals the ongoing star formation 
in the NH3 core and NW region. 



3.2 Distribution of the young population 

In the entire field, a total of 480 point sources are detected in at least 4 bands (including 2MASS JHKs 
bands). Using the combined photometric data, we construct the SEDs of the sources, which can help us 
to pick up local members from field stars. To do this, we classify the SEDs into 6 categories: 1. SEDs 



2 The online fitting tools can be accessed at http://caravan.astro.wisc.edu/protostars/ 
2 The WISE catalogue can be found at http://irsa.ipac.caltech.edu/Missions/wise.html 
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Fig. 3 Best fitting models for the two YSOs, Yl (left panel) and Y2 (right panel). The filled 
circles show the input fluxes. The black line shows the best fits. The dashed line shows the 
stellar photosphere corresponding to the central source of the best fitting model, as it would 
look in the absence of circumstellar dust (but including interstellar extinction). 



increasing monotonously towards longer wavelengths; 2. those having double peaks; 3. those with single 
peak at Ks or longer wavelengths; 4. those with single peak at H; 5. those decreasing monotonously 
towards longer wavelengths; 6. irregular SEDs. Fig.|4]shows some example SEDs of the above categories 
except for the type 6. 

According to the theory of star formation, the categories would represent an evolutionary sequence 
from 1 to 5. Even considering the extinction drops from the J-band through ch4, Type 1 SEDs (66/480) 
should peak at chl or longer wavelengths, so they might be Class I YSOs or still earlier; Type 2 SED 
sources (35/480) could be Class I or II objects with the presence of disks. Type 3 SEDs (108/480) are 
probably T Tauri objects with weak disk emissions; Type 4 SEDs (135/480) are more evolved, but could 
be still YSOs with infrared color excesses; Type 5 SEDs (35/480) may be local YSOs or field stars; For 
the last category, the fluxes measured are generally small with large error-bars. They might be YSOs or 
may be some line-emitting knots such as H2 or PAH features, or field stars or even false-detections. We 
cannot infer their properties so we ignore this kind in the following analysis. 
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Fig. 4 Typical SEDs of the 5 categories. The numbers correspond to our internal photometric 
identifications. The source No. 362, 509, 513, 563 and 1339 show typical SED type 1-5, 
respectively. We note however, not all sources are detected in all bands as the examples shown 
here. 



Fig. [5] shows the color-color diagram of the point sources detected in all four channels. It is not 
surprising to see that Type 1 sources are located in the Class 0/1 block or in the upper-right corner of the 
Class II block. Type 2 sources are in the Class II block while Type 3 sources are in the Class II and III 
block. Type 4 and 5 sources are mostly in the Class III block. Such a coincidence indicates the scheme 
of evolutionary sequence we proposed above is somewhat equivalent to the Class I— III sequence. We 
note here, it is difficult to tell the virtue of one scheme over the other. Our proposed scheme contains 
the information between 1 and 10 /im, while the color-color plot gives more accurate color index in a 
smaller wavelength range. 

To illustrate the distribution of the YSOs in the region, we make a pseudo-color image as follows: 
We use red points to denote the location where Type 1 SED YSOs are; same procedures are done for 
other types while Type 2 and 3 YSOs are marked by yellow points, Type 4 by green points, and Type 5 
by blue points. Such an image is then convolved by a 45"x45" kernel. The resultant image is presented 
in Fig. 

As expected, the red colors representing the youngest population that have been detected yet are 
located in the NH3 core, stretching southward along the edge of the shell-like structure surrounding the 
NIR cluster. Another red pea k is located in the NW r egion where two H2 jets are found. In Fig. [6] we 
overlay the 850 /im contours ( Di Francesco et al.|2008| 3 on the YSO distribution image. Interestingly, 
the red peaks coincide with the 850 /Ltm continuum peaks ( |Klein et al.|2005[ ) very well. This shows that 
the mm continuum emission is associated with the very young population there. This result confirms 
our previous conclusion that a cluster with two subgroups is just forming. Taking the NIR cluster into 
consideration, the whole region shows a star forming sequence that the NIR cluster forms first, and the 
expanding shell surrounding the NIR cluster may triggered another round of star formation in the NH3 
and NW areas. 

Using the JCMT archive data, which were obtained by Di Francesco et al. (2008 1, and assuming 
optically thin dust emission, a dust temperature of 20K, gas-to-dust ratio of 100, grain size of 0.1 /im, 
grain mass density of 3 g cm -3 , and a grain emissivity index of 2 (corresponding to n s» 0.3 for compar- 



3 The 850 fj,m image can be downloaded at 



http://wwwl.cadc-ccda.hia-iha.nrc-cnrc.gc.ca/community/scubalegacy/ 
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ison o f|Ossenkopf & Henning|1994| l, we adopt the equations outlined in Hildebrand]( |1983| l and |Beuther| 



et al. (2005 1 to estimate the gas masses. The final values are ~ 1100 M Q and ~ 150 M Q in the NH 3 



core and NW region, respectively. Meanwhile, the masses of the stellar contents are estimated by SED 
fittings to be ~ 240 M and ~ 50 M , respectively. The overall star formation efficiencies are then ~ 
18% and ~ 25% in these two subregions. 



4 CONCLUSION 



We have carried out a deep observation of the star forming region AFGL 5157 using the Spitzer/IRAC. 
PAH emissions are detected to be extended throughout the field. No destruction of PAH emissions 
suggests that there is no very high mass star in the region. The images in the four channels reveal a large 
number of point sources that are not detected in the NIR. By analyzing the SEDs of the point sources in 
the NIR and IRAC four channels, we have identified a deeply embedded cluster, which is split into two 
subgroups. These two subgroups contain a number of point sources with SEDs increasing toward longer 
wavelengths. They congregate toward the 850 fim continuum peaks (Di Francesco et al. 2008; Klein et 



al.|2005 I. The spatial coincidence between the SED I sources and the sub-mm peaks suggest that a new 



cluster is just forming there. The star formation efficiencies in the two subregions are estimated to be 
18% and ~ 25%, respectively. 

In the two subregions we detected a large number of objects with enhanced emission in ch2. Some 
of them have, but some have no counterparts in the NIR H2 narrow band image ( Chen et al.|2003[ ). We 
have identified several bipolar jets in the NH3 core region; two of them are associated with CO bipolar 
outflows. The SED fittings and jet kinematics both suggest their driving sources are very young. One of 
the driving sources is not detected even in the IRAC images, but shows faint mm continuum emissions 
(Font anTet al.||2009] >. Probably it has just started the forming process. Two jets that are presumably 
originated from the NW region are detected. The widespread H2 knots also support our conclusion that 
a cluster is forming recently. 
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Fig. 5 The [5.8]-[8.0] vs [3.6]-[4.6] color-color diagram of the point sources detected with 
signal-to-noise ratio better than 5 (photometric error less than 0.2 mag). The boundaries of 
Class 0/1 - III regimes is adopted from |Allen et al. | (|2004). 
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Fig. 6 Pseudo-color image representing the distribution of YSOs in the field. The red color 
represents for the Type 1 SED objects; the yellow for the Type 2 and 3; the green for the Type 
4 and the blue for the Type 5 (see text). The solid contours are SCUBA 850 /im continuum 
image ( Pi Francesco et al.|2008 Klein et al.|2005 1. 



